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Abstract: The paper reports on the coupling of Parity-Time (PT)- 

symmetric whispering gallery resonators with realistic material and 
gain/loss models. Response of the PT system is analyzed for the case of low 
and high material and gain dispersion, and also for two practical scenarios 
when the pump frequency is not aligned with the resonant frequency of 
the desired whispering gallery mode and when there is imbalance in the 
gain/loss profile. The results show that the presence of dispersion and 
frequency misalignment causes skewness in frequency bifurcation and 
significant reduction of the PT breaking point, respectively. Finally, as 
coupled WGM resonators are inherently lossy structures, we show that 
unbalancing the gain/loss in resonators is required to compensate for 
inherent loss of the structure and achieve improved PT properties. 
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1. Introduction 

Photonics is emerging as a popular practical platform for the exploration of Parity-Time (PT)- 
symmetric systems characterized by balanced loss and gain and having a threshold point at 
which real eigenfrequencies of the system coalesce to become complex conjugates IM.The 
existence of this threshold point is essential to the unique properties of PT-structures, such as 
unidirectional invisibility and simultaneous lasing and absorption i2im . This phenomena opens 
new avenues for the realization of practical devices such as lasers, optical memory, optical 
switches and logic-gates GHm . To date, PT-symmetric structures based on Bragg gratings and 



coupled optical systems have been investigated both theoretically ||2H4l[8l fT2l - [T9l and exper¬ 
imentally ll20lB5l . Recently, a PT-symmetric system based on two coupled microresonators 
and two fiber-taper waveguides has been experimentally demonstrated and shown to exhibit 
direction-dependent behavior at a record low power of 1/rW 112111251 . This is primarily at¬ 
tributed to strong field localization and build up of energy in the resonant whispering gallery 
modes M211125 1 . and has further strengthened the argument for using resonant structures rather 
than waveguides as building blocks of PT-symmetric systems. In contrast to PT symmetric cou¬ 
pled waveguide systems where the eigenmodes are purely real below the threshold point, the PT 
symmetric coupled microresonators have complex eigenfrequencies below the threshold point 
due to inherent radiation losses EHiniini. 

In this paper we investigate the fundamental properties of the PT resonant system based 
on two coupled whispering gallery resonators within the context of both realistic material 
properties and practical operating constraints. In particular we discuss how practical disper¬ 
sive properties of material gain and loss that satisfy the Kramers-Kronig relationship affect the 
performance of microcavity-based PT resonant structures. Our surprising conclusion is that ac¬ 
counting for large, yet realistic, levels of dispersion preserves the essential threshold-behaviour 
predicted by completely PT-symmetric dispersionless models, while more moderate levels of 
dispersion can completely change the character of the response of the system to increasing gain 
and loss. In particular, when there is moderate dispersion the gain and loss materials respond 
differently to frequency shifts in such a way that sharp threshold points give way to gradual 
changes over a range of parameters. When dispersion is increased further, the response reverts 
to threshold behaviour of the type seen in non-dispersive PT-symmetric systems, albeit with 
some breaking of detailed quantitative symmetry. 

Our recent work on dispersive PT-Bragg gratings has shown that material dispersion limits 
the unidirectional invisibility to a single frequency which is in stark contrast to previous results 
that assumed idealized gain/loss profile in order to demonstrate wideband unidirectional behav¬ 
ior M191I26M . In this paper, the performance of the microresonator-based PT system is analyzed 
for practical scenarios involving: a) frequency mismatch between the cavity resonant frequency 
and the gain pump frequency and, b) imperfect balance of the gain and loss in the system. The 
analysis of the microresonator-based PT system is achieved using an exact representation of 
the problem based on boundary integral equations and explicit analytical results are given for a 
weakly coupled system using perturbation analysis Et). We concentrate on the weakly-coupled 
limit in our detailed calculations because that captures the essential properties of the threshold 
behavior of the PT-symmetry while allowing simple analytical calculations to be used. Finally, 
real-time field evolution in a two microresonator PT-symmetric system is analyzed for differ¬ 
ent levels of dispersion using the numerical time domain Transmission Line Modelling (TLM) 
method 

2. PT symmetric coupled microresonators 

In this section we describe the theoretical background of a PT-symmetric system based on two 
coupled microresonators. The system, in which both microresonators have radius a and are sep¬ 
arated by a gap g, is illustrated schematically in Fig. [T] The active and passive microresonators 
have complex refractive indices hq and iii respectively, that are typically chosen to satisfy the 
PT condition hq = n£, where * denotes complex conjugate, n= {n' + jn"), and n' and n" rep¬ 
resent the real and imaginary parts of the refractive index. In practice, localized gain might 
be achieved by means of erbium doping and optical pumping of the active microresonator, 
while masking the lossy microresonator as in Il20lj22ll25l . Both resonators are assumed to be 
surrounded by air. 

The refractive index of dispersive materials is frequency dependent but must also satisfy the 



Fig. 1. Schematic of two coupled cylindrical microresonators or radius a and separated by a 
distance g. Microresonators with gain and loss are denoted by fxRc and flRi, respectively. 


causality property between the real and imaginary parts of the material refractive index I30II31I . 
The material properties are conveniently modelled by assuming a dielectric constant that uses 
a Lorentzian model for dispersion as in ll32l 

s(o)) = e ■ ^0 ^ I 1 

Sr (0) Soc J2£i^c0\l+j{C0 + (Oa)T: I + jicO - C0a)T: 

Here Soo denotes the permittivity at infinity, cOa denotes the atomic transitional angular fre¬ 
quency, T is the dipole relaxation time and Oo is related to the conductivity peak value that is 
set by the pumping level at cOa- The time-varying component has been assumed to be of the 
form e-'®' and therefore Ob > 0 denotes loss while Ob < 0 denotes gain. The parameter T controls 
the degree of dispersion, with T = 0 corresponding to a dispersion-less system. Throughout this 
paper, the frequency-dependent refractive index is expressed as n = Sr{co) and the material 
gain/loss parameter is expressed using the imaginary part of refractive index as 7 = (On . 



3. Analysis of inter-resonator coupling in the frequency domain 

We now provide an analysis of coupling between resonators based on boundary integral meth¬ 
ods. This approach is particularly suited to perturbative approximation of the coupling strength 
in the weak coupling limit but also provides an efficient platform for exact calculation when 
coupling is strong. The calculation is based on an approach used in lIZTll to describe coupling 
between fully bound states in coupled resonators and optical fibers, but adapted here to allow 
for radiation losses. It is also similar to methods used in ll33l435]l . 


3.1. Notation and assumptions 


We assume resonators of radius a, uniform refractive index and TM boundary conditions. Then 
the mode taking the form Xj/i — {Jm{nikr)/Jm{nLka))e^'^^ inside the isolated lossy resonator is 
such that the solution and its normal derivative on the boundary of the resonator can be written 
as 


dWL 

= 

on 


F^Wl, 


( 2 ) 


where 


pL _ (^) 

“ Jm{z) 


and z = n^ka, 


(3) 


where, k is the free-space wave number and y/Q and being dehned similarly for the gain 
resonator. The treatment of coupling in the remainder of this section can be used for other 
circularly-symmetric resonators such as those with graded refractive index or with different 
boundary conditions, as long as an appropriately modihed is substituted in (| 2 |i. 










3.2. Exact solution using boundary-integral representation 


An exact boundary integral representation of the coupled problem is conveniently achieved 
by applying Green’s identities to a region fi which excludes the resonators, along with an 
infinitesimally small layer surrounding them (so that the boundaries Bq and Bl of the resonators 
themselves lie just outside H). In we assume that the refractive index takes the value no = l, 
so that the free-space Green’s function is 


Go(x,x') = -^//oWx-x'|), (4) 

where Ho{z) = Jo{z) — jYo{z) denotes the Hankel function of the second kind (and the solution 
is assumed to have time dependence e^®')- Then, applying Green’s identities to the region fi 
and assuming radiating boundary conditions at infinity leads to the equation 


0 = 



dGojxy) 

dn' 



(5) 


when X lies on either B^ or Bq (and therefore just outside of f2). 

We now expand the solution on the respective resonator boundaries as Fourier series, 

re = 1 ^ 0 !,and V4 = ^ 

m m 


in the polar angles 9q and 9l centered respectively on the gain and lossy resonators, running 
in opposite senses in each resonator and zeroed on the line joining the two centers. The corre¬ 
sponding normal derivatives can be written 


dWa 

dn 


1 m ^ 


and 


dWG 

dn 


- ;m 0 G 


(7) 


Using Graf’s theorem ll^ to expand the Green’s function Go(x,x') analogously in polar coor¬ 
dinates on each boundary, the integral equation (|5]l, evaluated separately for x on Bq and on Bq, 
leads to a pair of matrix equations 


+ = 0 


( 8 ) 


Here, 


«« = 

1 

(yp 

r^.G 

and y = 



^m+1 


G^m+1 


(9) 


are Fourier representations of the solution on the boundaries of the gain and lossy resonators 
respectively. The matrices and are diagonal with entries 





timijY) ) 


where u = ka, 


( 10 ) 


and provide the solutions of the isolated resonators. The matrices and describe cou¬ 
pling between the resonators. The matrix has entries of the form 


/^GL 

^Im 


Jl{u)Hi+m{w)Jm{u) 



ULJm{u) \ 
dmi^) ) 


( 11 ) 


where u = ka, w = kb and b is the center-center distance between the gain and lossy resonators. 
The matrix is defined analogously by exchanging the labels G and L. 












3.3. PT-symmetry in the exact solution 

The system ® can be presented more symmetrically by using the scaled Fourier coefficients 


f/L _ j / N (pL _ uJ'm{u) \ , 

«m - Jm{u) \F„ j 

(along with an analogous definition of a,^). Then ® can be rewritten 

b^a^ + Ca^ = 0 

CaP^b^a^ = 0 , 


where the diagonal matrices have entries 


b^P = - 7 
^mm J 




Jm{u)Fp’^ - uJ'^{u) 


where u = ka, 


and the same (symmetric) matrix C, with entries 




( 12 ) 


(13) 


(14) 


(15) 


couples solutions in both directions. 

We have included an overall factor of —j in these equations to emphasise an approximate PT- 
symmetry that occurs when no = n*^. Then, in the limit of high-Q whispering gallery resonances 
for which we may approximate 

jH,n{u) ~Ym{u) and jHi+^{u) ~Yi+„^{u), (16) 

the matrices in (fOT l satisfy the conditions 

{b^y^b^ and C*~C (17) 

which are a manifestation of PT symmetry of the system as a whole; deviation from these 
conditions is due to radiation losses. 

3.4. Perturbative weak-coupling approximation 

The system of equations ( fTST l can be used as the basis of an efficient numerical method for 
determining the resonances of the coupled system with arbitrary accuracy. In practice, once the 
gap g = b — 2a between the resonators is wavelength-sized or larger, a truncation of the full 
system to a relatively small number of modes suffices to describe the full solution. 

In the limit of very weak coupling we may achieve an effective perturbative approximation 
by restricting our calculation to a single mode in each resonator. We consider in particular the 
case of near left-right symmetry in which 


nG~«L- (18) 

PT symmetry is achieved by further imposing ng = but for now we allow for the effects of 
dispersion by not assuming that this is the case. We build the full solution around modes for 
which 

\I/±^\I/g±Wl, (19) 

where y/c und xj/i are the solutions of the isolated resonators described at the beginning of this 
section. We use a single value of m for both y/e and xj/i and in particular approximate the global 




mode using a chiral state in which the wave circulates in opposite senses in each resonator. That 
is, we neglect the coupling between m and —m that occurs in the exact solution. 

Then a simple perturbative approximation is achieved by truncating the full system of equa¬ 
tions (foi l to the 2x2 system 

M ( “r” 1=0, where M=( 1 . (20) 

\ OCmm / \ Cmm / 

Resonant frequencies of the coupled problem are then realised when 

0 = det A/ = 1211 

In the general, dispersive and non-PT-symmetric, case this reduces the calculation to a semi- 
analytic solution in which we search for the (complex) roots of the known 2x2 determinant 
above, in which matrix elements depend on frequency through both k = co/c and n = n{(o). 

3.5. Further analytic development of the perturbative solution 
To develop a perturbative expansion we let 

rfi = — (fP f)^ ) and — — (f)^ — fP ) 1221 

^mm 2 aiiu ‘-'mm y-^mm ‘-^mmj 

(and note that in the high-Q-factor PT-symmetric case, ~ (5^)*, both and are 
approximately real). We assume that both and Cmm are small and comparable in magnitude, 
and expand angular frequency 

Aoto 

coi .2 = (Oo±-^ + --- (23) 

about a real resonant angular frequency of an averaged isolated resonator satisfying 

Dim{(Oo) = 0. (24) 


Then to hrst order the coupled resonance condition becomes 

0 = detM = {(Oof + - Cmm{(0of + ■■■ 

from which the angular frequency shifts can be written 

Aft)o ^ \/Cmm (^ ^rnm ( ^ 

DOjicoo) ’ 

where denotes a derivative of Dmm{<^) with respect to frequency. 

We then arrive at a simple condition 


(25) 


(26) 


CnmiCOo)^ = D^,n,niO)of 

for threshold (at which Acuo and the two resonant frequencies of the coupled system collide). In 
the PT-symmetric case, where Cmm and Dmm are approximately real (and whose small imaginary 
parts represent corrections due to radiation losses), we therefore have a prediction for a real 
threshold frequency. 




4. Results and discussions 


In this section, the impact of dispersion on the performance of the PT coupled microresonators 
is analyzed. Frequency mismatch between the resonant frequency of the microresonator and 
gain pump frequency is investigated for practical levels of dispersion and the practical impli¬ 
cations of a slight unbalance between the gain and loss in the system are investigated. We 
conclude the section with an investigation of how coupling between resonators manifests itself 
in the time development of solutions. 

4.1. Ejfects of dispersion on threshold behavior in the frequency domain. 

For all cases, weakly coupled microresonators are considered, the coupled resonators with 
a dielectric constant e«, = 3.5 ll^ have radius a = 0.54/rm and are separated by distance 
g = 0.24/rm. Transverse-magnetic (TM) polarization is considered and operation at two dif¬ 
ferent whispering-gallery modes is analysed, namely a low Q-f^ctor mode (7,2) and a high 
Q-factor mode (10,1). The corresponding isolated resonator resonant frequencies are respec¬ 
tively = 341.59THZ and = 336.85THz, with g-factors = 2.73 x 10^ and 

g(i04) = 1.05 X lOl 

Figure |2] shows the real and imaginary part of the eigenfrequencies Oi and ©2 of the PT- 
symmetric coupled microresonators with balanced gain and loss, 7o = —7g = Yl, and is depicted 
as a function of the gain/loss parameter 70 = \cOon (©o)| for both the low and high Q-factor 
modes. The gain and loss are assumed to be tuned to the resonant frequency of an isolated 
microresonator, i.e. cOa = COo = 2;r/o. Three different levels of dispersion, controlled by the 
parameter T defined in Sec. 2, are considered. These are cOat: = 0 corresponding to the case 
of no dispersion, cOaT = 212 taken from lf3^ to exemplify the case of high dispersion and 
COaT = 0.7 to exemplify the case of low dispersion. 

Figure I2a,b) shows the frequency splitting of the real and imaginary part of the complex 
eigenfrequencies for the case of no dispersion. In the absence of gain/loss, where 70 = 0, the 
supermodes beat at a rate corresponding to the frequency differences ©1 — ©2 = 3.823 rad/ps 
and 1.164 rad/ps for the (7,2) and (10,1) modes respectively. Figure|2a) indicates that operation 
in a higher Q-factor mode results in weaker coupling between the microresonators compared to 
the case of operation in the lower Q-factor mode. Increasing the gain and loss in the system, de¬ 
creases the beating rate and the supermodes coalesce at the threshold points of 70 = 6.86 rad/ps 
and 2.1 rad/ps for the low and high Q-factor modes of operation respectively, confirming that 
the high -2 factor mode has a lower threshold point II 2 TI . In the case of operation in the low 
2-factor mode, the eigenfrequencies shown in Fig. |2b) have a significant constant and posi¬ 
tive imaginary part before the threshold point, which is a consequence of the higher intrinsic 
losses due to radiation in that case. The corresponding imaginary part is insignificant in the 
case of the high Q-factor mode, for which radiation losses are much smaller. Furthermore it is 
noted here that the coupled system first starts to lase, i.e. one of the eigenfrequencies satisfies 
Im(©i 2 — ©b) < 0 , only when operated significantly beyond the threshold Jq — I rad/ps for the 
low 2 -factor operation while this onset occurs immediately after the threshold point in the high 
2 -factor case. 

Figure|2c,d) shows the real and imaginary parts of the eigenfrequencies for the case of strong 
dispersion, corresponding to the parameter values (Oa'C — 2\2 taken from lf32l . These are again 
shown for both high and low 2-factor modes. It is noted that the threshold point for the low 
2-factor mode is reduced from 70 = 6.86 rad/ps to 70 = 6.47 rad/ps in this case while for the 
high 2 -factor mode it remains unchanged at 2.1 rad/ps (compared to the case of no dispersion). 
Below the threshold point the imaginary parts of the eigenfrequencies are not constant, but are 
instead skewed towards a lossy state with positive and increasing imaginary part. Extension 



7o = Won”(wo) / rad ps" 



70 = uian{u>a) / rad ps" 


Fig. 2. Frequency bifurcation of PT-coupled microresonator with a balanced gain (yg = 
—yo) and loss (y/, = yo) as a function of gain/loss parameter at the peak of pumping beam 
yo = (Oafi (cOa) for three different dispersion parameters, (a,b) cor = 0, (c,d) cor = 212 and 
(e,f) cor = 0.7 


beyond the threshold point shows that the imaginary parts of the eigenfrequencies do not split 
evenly and are also skewed towards overall loss, implying that in the highly dispersive case the 
eigenfrequencies both are complex but no longer complex conjugates after the threshold point. 

The real and imaginary parts of the eigenfrequencies for the case of low levels of dispersion, 
for which we take cOaT = 0.7, are shown in Fig.|2je,f). Figure|2|e) shows that there is no clear 
threshold point in this case: the imaginary parts split for very low value of the gain/loss param¬ 
eter 7o, with no sharp point of onset. The appearance of a threshold point typically associated 
with PT-behavior is lost and the eigenfrequencies are always complex valued. 

The key conclusion to be made from Fig. |2] is therefore that PT-like threshold behavior is 
observed in the cases of no dispersion and of high dispersion, but not for cases of intermedi¬ 
ate dispersion. While there is some skewness in the high-dispersion case, which amounts to a 
quantitative deviation from strict PT-symmetry, there is an essential qualitative similarity to the 
dispersionless case in which there appears to be a sharp threshold. By contrast, in the case of 
intermediate dispersion there is no sharp transition point and the imaginary parts of the two 
frequencies begin to diverge from the beginning. 

To further investigate and explain this phenomenon, we examine the dependence of the real 




























Fig. 3. (a) Impact of dispersion to the real part of material at atomic transitional angular 
frequency cOa due to the presence of gain and loss for different dispersion parameters; (b) 
Contrast between the real part of eigenfrequencies of PT-coupled microresonators for two 
different gain/loss parameter, i.e. 70 = 7.5 rad/ps for (7,2) and 2.54 rad/ps for the (10,1) 
mode as function of dispersion parameter T. 



Fig. 4. Frequency bifurcation of coupled microresonators with balanced gain and loss as 
function of gain/loss parameters ya, for two different atomic transitional frequencies cOa = 
iTtifo + 5) with 5 = -0.1 and 0.1 THz. 


part of the complex refractive index on the dispersion parameter Ico^T. This dependence is plot¬ 
ted in Fig. [3a) for the cases of both gain and loss, for which we respectively take ffq = ±2eo(bc7 
and (Oa — (Oo- Figure|3a) shows that the real parts of the refractive indices behave differently 
for the cases of loss and gain in the system, with the maximum difference occurring when 
T = 1/ (2a)a)- However, in two limiting cases T = 0 (dispersion-less system) and T —?> 0 ° (strong 
dispersion), the real parts of the refractive index converge. This means that the PT condition 
hq = can only be satisfied accurately for the cases of no dispersion and of high dispersion. 
For the case of intermediate dispersion there is necessarily some discrepancy between the real 
parts of the refractive indices of the resonators. 

Figure |3b) shows the minimum difference in the real parts of the two eigenfrequencies 

















for different dispersion levels and operated at a fixed value of the gain/loss parameter, i.e. at 
7o = 7.5 rad/ps for the low Q-factor and at 2.54 rad/ps for high Q-factor modes of operation. 
These values of the gain/loss parameter are chosen to lie above the expected threshold so that 
qualitatively PT-like behavior would imply eigenfrequencies with a common real part. Figures 
Oa) and |3b) confirm that the maximum difference between real parts of the two refractive 
indices coincides with the maximum deviation from PT-like threshold behavior, where the dif¬ 
ference between real parts of the eigenfrequencies is greatest. This result further confirms the 
fact that realistic levels of dispersion preserve the essential features of PT behavior. 

Having conhrmed that realistic levels of dispersion preserve PT behavior, Figure|4]considers 
a practical scenario in which there is high dispersion (OaT — 212 and a frequency mismatch be¬ 
tween the resonant frequency and the gain/loss atomic angular frequency. The material atomic 
frequency is defined to be (Oa — 27r(fo + S), where d is the mismatch parameter. The structure 
is operated with balanced gain and loss, i.e. Yg = —Yl and two values are assumed for the fre¬ 
quency mismatch, namely 5 = —0.1 and 0.1 THz. Figure |4|a,b) shows the results for the low 
2-factor mode (7,2) and Fig. |4jc,d) for the high-Q factor mode (10,1). In both cases there is 
no sharp threshold point for the real parts of eigenfrequencies and the imaginary parts begin 
to diverge at low gain/loss values. Neither are the imaginary parts symmetrically placed about 
a branching value. This result conhrms the fact that PT behavior is preserved only when the 
angular transitional frequency of the dispersive gain/loss prohle is aligned with the resonant 
frequency of the microresonators. If that is not the case, the frequency misalignment causes the 
coupled system to continue to beat after a threshold region. 




Fig. 5. Complex eigenfrequency in a PT-coupled microresonator system with variable gain 
and fixed loss shown as a function of gain parameter \Yg\< dispersion parameter (no-r = 212 
and shown for 3 different fixed loss value, i.e. = 5.565, 6.4281, and 7.291 rad/ps. 

Another practical scenario is considered in Fig. |5] where the gain and loss are not balanced, 
i.e. iiRl has a loss Yl while iiRg has a gain Yg- Figure|3a,b) shows the real and imaginary parts 
of the eigenfrequency for three different values of loss namely, Yl — 5.565, 6.4281 and 7.291 
rad/ps which correspond to values below, at, and above the threshold point of a PT symmetric 
structure with balanced gain and loss respectively. The low g-factor mode is considered with a 
practical dispersion parameter of OcrT = 212 as taken from ll32l . Interestingly, we now observe 
that the PT threshold point can also exist for structures with unbalanced gain/loss as shown 
by the plots for Yl — 5.565 rad/ps and Yl — 7.2910 rad/ps in Fig. |5] In the former case, the 
PT threshold is increased and in the latter case the PT threshold is decreased compared to the 
PT threshold of the balance structure. Of special interest is the observation that increasing loss 
results in the reduction of the PT threshold which consequently reduces the levels of gain at 
which lasing occurs. This counter-intuitive principle of switching lasing on by increasing loss 
has been experimentally demonstrated in ll22l where a metal probe is used to enhance loss in 
the lossy microresonator. 
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Fig. 6. (a) Spatial electric field distribution of the coupled microresonators operated (7,2) 
mode. The black dashed line denotes the monitor line. The temporal evolution (b) and 
spectra (c) of the field on the monitor line are shown in the absence of gain and loss. 


4.2. Real time operation ofPT symmetric coupled microresonators 

In this section the real-time operation of the PT-symmetric coupled microresonators is demon¬ 
strated for different levels of dispersion. For this purpose, the two-dimensional (2D) time- 
domain Transmission Line Modeling (TLM) numerical method is used. A more detailed de¬ 
scription of the TLM method is given in II28II29I and the implementation of general dispersive 
materials for PT-symmetric Bragg gratings is demonstrated in ifTolIT^ . In each of the simu¬ 
lations shown in this section, the low Q-factor (7,2) mode is excited by a very narrow-band 
Gaussian dipole located in pRq whose frequency is matched to the resonant frequency of this 
mode. Depending on the levels of gain and loss, and their relation to the threshold points, we 
find in practice, however, that small unintentional initial excitations of the high Q-factor ( 10 , 1 ) 
mode may grow to become a significant feature and even dominate the evolved state. In all 
cases we find that the TLM simulations are consistent with the frequency-domain calculations 
provided in the previous section and in fact have been used to independently validate the pertur¬ 
bation analysis results presented in Figs.|2ll5] although the detailed calculations are not reported 
here. 

We begin with the case of the evolution from the low Q-f^ctor mode using a model with no 
dispersion. Figure | 6 t a) shows the spatial electric field distribution of coupled microresonators 
with no gain and loss (70 = 0) and operating at the resonant frequency of the low Q-factor (7,2) 
mode. The black dashed line denotes the monitor line on which the electric field is observed 
during the TLM simulation. Parts (b,c) of Fig.| 6 ]show the temporal evolution and the spectra of 
the electric field observed along the monitor line for the case of no gain and loss. 

The case of no gain and loss, reported in Fig. 6 (b), shows a typical oscillation of the electric 
field between the microresonators having a regular beating pattern in which maximum intensity 
being observed in one microresonator corresponds to minimum intensity being observed in the 
other. Figure 6 (c) shows the frequency content of the modes, indicating the presence of two 
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Fig. 7. The temporal evolution and spectra of the field on the monitor line are shown for 
different gain/loss parameters, (a,b) for yo = 4-3 rad/ps and (c,d) for Yq = 7.5 rad/ps with a 
negligible dispersion parameter using the TLM method. 


(1 2 ) 

resonating frequencies centered around /q ’ , in agreement with Fig. 2(a). 

The real-time performance of PT-coupled microresonators with balanced gain and loss Yg — 
~7o^ Yl = 7o with no dispersion is depicted in Fig.|2l The temporal evolution and the spectra of 
the electric field are observed along the monitor line for two levels of gain and loss: levels of 
gain and loss 70 = 4.3 rad/ps, i.e. lower than the threshold point of the low Q-f^ctor mode but 
beyond the threshold of the high Q-factor mode in Fig.|7](a and b) and levels of gain and loss 
7o = 7.5 rad/ps, i.e. lying above the threshold points of both modes in Fig.|2](c and d). 

Figure [TJa) shows the energy beating between the microresonators with gain/loss 70 = 
4.3 rad/ps, set below the threshold point of the low-Q-factor mode. It is noticeable that beating 
between microresonators is reduced and is no longer periodic. Additionally, Fig.|7Ja) indicates 
the presence of additional modes, observable at later times. Frequency analysis of the fields 
is given in Fig. |2tb) and shows an additional peak at 336.85 THz, which corresponds to the 
resonant frequency of the mode (10,1), explaining the high frequency beating in Fig.|7ja). Re¬ 
ferring to Fig. |2ta) it can be seen that at 70 = 4.3 rad/ps, the (10,1) mode is operating above 
its threshold point and thus experiencing amplification, whilst the (7,2) mode is still below its 
threshold point. 

A further increase in gain/loss to 70 = 7.5 rad/ps, in Fig.ITJc), shows an exponentially grow¬ 
ing field in the gain microresonator with no beating between the resonators and a stronger pres¬ 
ence of the high Q-factor mode as shown in Fig. [TJd). Referring to Fig.|2ja,b) it is confirmed 
that for 7 o =7.5 rad/ps both low and high Q-factor modes are operating above the threshold. 

Figure m shows corresponding results for the case of high dispersion with cUcrT = 212. We 
once again choose gain/loss parameters 70 = 4.3 rad/ps and 70 = 7.5 rad/ps but omit the case of 
zero gain/loss here. Figure [ 8 l a) shows a decaying beating pattern. The corresponding spectral 
analysis in fig.Etb) shows that the beating may be attributed to low Q-factor modal frequencies, 
indicating that the highly dispersed gain/loss profile has stabilized the operation of PT-coupled 
resonators system at a desired mode of operation. For operation with gain/loss parameter 70 = 
7.5 rad/ps, the temporal response in Fig.[ 8 jc) indicates an exponentially growing field with no 
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Fig. 8. Temporal and spectra of electric field along the monitor line for coupled PT mi¬ 
croresonators with balanced gain and loss parameters operated for (7,2) mode with prac¬ 
tical dispersion parameters cOa'^ = 212 t32\ and for two different gain/loss parameter, i.e. 
(a,b) 7o = 4.3 rad/ps and (c,d) 70 = 7.5 rad/ps. 


presence of high order modes. The spectrum in Fig. Od) shows a single peak at the resonant 
frequency of the (7,2) mode, confirming that the resonators are operating above the PT threshold 
point. The strongly dispersive gain/loss profile limits operation of PT-coupled microresonator 
system to the low Q-factor (7,2) mode only in this case. This again confirms the result that when 
the material atomic frequency is chosen to be at a desired resonant frequency, the PT-symmetry 
is limited to that particular mode only, as in the case of periodic potentials I19II31L 

Figure |9] investigates the real-time operation of the PT structure with unbalanced gain and 
loss shown in Fig. |5] Here, we apply the same scenario as in lf22l where the gain in the active 
microresonator is fixed at Yg = —7.053 rad/ps and the loss is varied in the passive microres¬ 
onator. The low Q-factor mode with practical dispersion parameters of cOaT = 212 is consid- 
erd ll^ . Figure |9ta) shows the electric field observed along the monitor line for the case of 
loss Yl = 5.565 rad/ps, i.e. more gain than loss in the system. It can be seen that there is a 
non-periodic and long temporal beating pattern. The field is not growing which indicates that 
the system is stable and is not lasing. The spectral decomposition shown in Fig. |9jb) is unable 
to distinguish the splitting of the resonant frequency is due to the limited spectral resolution of 
the Fourier transformation of the time domain simulation result. 

Figure |9jc) shows the temporal evolution of the field for the case of Yl = 6.4281 rad/ps. It 
can be seen that there is no beating between the resonators and no growing field, suggesting 
that structure is operating above the PT threshold point but before the lasing point which occurs 

Yg = —7.377 rad/ps (FiglSjb)). The spectral analysis shows only a single peak centered at 
AT.2) 

Jo 

Figure |9le) shows the temporal evolution for Yl = 7.291 rad/ps and Yg = —7.053 rad/ps, 
i.e. more loss than gain in the system. It can be observed that the field is growing with no 
beating between the microresonators, suggesting lasing action. This result is in agreement with 
observations in ll22ll where loss induced lasing is demonstrated. Corresponding spectrum is 


























Fig. 9. Temporal and spectra of electric field along the monitor line for coupled PT mi¬ 
croresonators with unbalanced gain and loss operated for (7,2) mode with practical disper¬ 
sion parameters cOat = 212 II32I . i.e. (a,b) Yl = 5.565 rad/ps, (c,d) Yl = 6.4281 rad/ps and 
(e,f) Yl = 7.291 rad/ps while the gain parameter is kept constant at Yg = -^7.053 rad/ps. 


depicted in Fig.|9tf) has a sharp peak centered at 

5. Summary and conclusion 

In the paper the impact of material dispersion on PT-symmetric coupled microresonators has 
been analyzed. It has been shown that the practical case of high dispersion preserves the require¬ 
ment for a PT structure. However our results shows that this is only the case when the material 
atomic frequency is aligned with the resonant frequency of the microresonator. This comes as a 
direct consequence of the Kramers-Kronig relationship which implies that changes in the imag¬ 
inary part of the refrective index caused the real part of the refractive index to change too. In ad¬ 
dition, we also demonstrate the principle of loss-induced lasing mechanism which is triggered 
by an early PT-symmetry breaking. Real-time operation of PT-coupled microresonators verifies 
that the dispersion due to the Kramers-Kronig relationship limits the operation of PT-coupled 
microresonators to a single frequency and hence forbids multi-mode PT-symmetry breaking. 



















































